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Abstract 

We propose a modified Balitsky-Fadin-Kuraev-Lipatov equation from the 
viewpoint of the resummation technique, which contains an intrinsic depen- 
dence on momentum transfer Q, and satisfies the unitarity bound. The idea 
is to relax the strong rapidity ordering and to restrict phase space for real 
gluon emissions in the evaluation of the BFKL kernel. It is shown that the 
power-law rise of the gluon distribution function with the small Bjorken vari- 
able x turns into a logarithmic rise, and that the predictions for the proton 
structure function F 2 (x,Q 2 ) are consistent with the HERA data. 
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1. Introduction 



It is known that the Balitsky-Fadin-Kuraev-Lipatov (BFKL) equation [|TJ 
sums leading logarithms ln(l/x), x being the Bjorken variable, which are 
produced from the reggeon ladder diagrams with the rung gluons obeying 
the strong rapidity ordering. This equation predicts a rise of the gluon dis- 
tribution function and thus a rise the proton structure function F2 involved 
in deep inelastic scattering (DIS) at small x, which have been confirmed by 
the recent HERA data 0. However, some controversies remain unsolved. 

Since the BFKL equation is independent of momentum transfer Q, its pre- 
dictions are insensitive to the variation of Q. On the contrary, the experimen- 
tal data for F 2 (x, Q 2 ) exhibit a stronger Q dependence ||: The rise at low x is 
slower for smaller Q, corresponding to so-called soft pomeron exchanges. The 
rise obtained from the BFKL equation is always larger, which corresponds 
to hard pomeron exchanges. To explain the data, the Dokshitzer-Gribov- 
Lipatov-Altarelli-Parisi (DGLAP) equation || , which sums single logarithms 
In Q produced from strong transverse momentum ordering, is then employed. 
The relevant splitting function P gg contains a term 1/x, which also gives a 
rise at small x. An alternative solution || is the Ciafaloni-Catani-Fiorani- 
Marchesini (CCFM) equation, which embodies the DGLAP and BFKL equa- 
tions ||. For this unified evolution equation, the Q dependence comes from 
the summation of the ladder diagrams with strong angular ordering. 

On the other hand, the increase of the gluon distribution function at 
x — > predicted by the BFKL equation is power-like, such that F2 and the 
DIS cross section a tot rise as a power of x, i.e., as a power of s for x = Q 2 /s, 
s being the total energy. This behavior does not satisfy the Froissart bound 
Ctot < const, x In s, and violates unitarity. Hence, the BFKL equation can 
not be the final theory for small x physics. Though it has been expected 
that the inclusions of next-to-leading ln(l/x) || and of higher-twist effects 
from the exchange of multiple pomerons |7j may soften the BFKL rise, the 
attempts have not yet led to a concrete conclusion. 

Recently, we have proposed a modified BFKL equation from the view- 
point of the resummation technique ||, which involves an intrinsic Q depen- 
dence 0. This Q dependence was introduced by cutting off the longitudinal 
component of the real gluon momentum at a scale of order Q, when evaluat- 
ing the BFKL kernel. The motivation is that it is unlikely for a real gluon to 
carry an infinite momentum. The gluon distribution function derived from 
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this Q-dependent evolution equation gives predictions of F2, which are in 
good agreement with the data. Unfortunately, the predictions still rise as a 
power of x. In [[TTJ we observed that relaxing the strong rapidity ordering 
for real gluon emissions results in a destructive correction, and the power- 
law rise of the gluon distribution function turns into a logarithmic rise. The 
unitarity bound is thus satisfied. 

In this letter we shall combine the above two modifications to derive a new 
small-x evolution equation, which depends on Q and satisfies the unitarity 
bound at the same time. We demonstrate that the predictions for F2 are 
consistent with the current data for x > 10~ 4 , and that the logarithmic rise 
of F2 turns on as x < 10~ 4 . 



2. The formalism 

We refer the detailed derivation of the the BFKL equation using the 



resummation technique to [|TTJ. Here we summarize only the basic idea. The 
unintegrated gluon distribution function F(x, kr), describing the probability 
of a parton carrying a longitudinal momentum fraction x and transverse 
momenta k^, is defined by 

F(x,hr) = —f—f d ^e-^v~-^y T ) 
K ' ' p+ J 2tt J 4tt 

xW(p,a\F+(y-,y T )F^(0)\p,a) , (1) 

Z a 

in the axial gauge n ■ A = 0, n being a gauge vector with n 2 7^ 0. The ket 
\p, a) denotes the incoming proton with the light-like momentum p^ = p + 5 lJ-+ 
and spin a. An average over color is understood. F+ is the field tensor. 

For a fixed parton momentum k + , F varies with p + implicitly through 
x = k + /p + , and we have —xdF/dx = p + dF/dp + . That is, the differentiation 
of F with respect to x is related to the differentiation with respect to p + . 
Because of the scale invariance in n of the gluon propagator —iN^ u (l)/l 2 with 

N K» = g v» + n 2^^ ; 2 

F must depend on p and n via the ratio (p ■ n) 2 /n 2 . Hence, there exists a 
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chain rule relating p + d/dp + to d/dn, 

d „ n 2 d „ 
dp + v ■ n drip 

V/3 = 5/3 + being a vector along p. The operator d/drip applies to a gluon 
propagator, giving 

d -N vv ' = l —U v N Pv> + r'N^) . (4) 



drip n ■ I 

The loop momentum l v (l v ') contracts with the vertex the differentiated gluon 
attaches, which is then replaced by a special vertex vp = n 2 vp/(v ■ nn ■ I). 
This special vertex can be read off from the combination of Eqs. (0) and @. 

Summing the diagrams with different differentiated gluons, and employ- 
ing the Ward identity for the contraction of l v il v ), those diagrams embed- 
ding the special vertices cancel by pairs, leaving the one in which the special 
vertex moves to the outer end of the parton line 0. We then obtain the 
derivative, 

- x -^-F(x,k T ) = 2F(x,k T ) , (5) 
dx 

where the new function F contains one special vertex [11]. The coefficient 



2 comes from the equality of the new functions with the special vertex on 
either of the two parton lines. 

Next we factorize the new function into the convolution of the subdiagram 
involving the special vertex with the original gluon distribution function. The 
resultant factorization formula will lead to the small-x evolution equation. 
The kinematic approximation, i.e., the rapidity ordering of radiative gluons, 
is specified only at the stage of computing the subdiagram, which is then 
identified as the evolution kernel. Therefore, in our approach it is easier to 
handle the approximation applied to the evolution kernel, and to investigate 
its effect on the behavior of the gluon distribution function. 

The only important region of the loop momentum / flowing through the 
special vertex is soft, in which the subdiagram containing the special vertex 
can be factorized. Including the color factor, the factorization formula is 
written as 0, [10 



F(x,k T ) = iN c g 2 J ^N^(l) V -^[2 m 6(l 2 )F(x + l+/p + ,\kT + lT\) 



(27T) 
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+ e ^~® F(x,hT) 



p 



(6) 



The first term in the brackets corresponds to the real gluon emission, where 
F(x + l + /p + , |k T + 1 T |) implies that the parton coming out of the proton 
carries the momentum components xp + + / + and + It in order to radiate 
a real gluon of momentum I. The second term corresponds to the virtual 
gluon emission, where the 6 function sets the upper bound of It to kr to 
ensure a soft momentum flow. 
The integration over l~ gives 



/■'(.r. h;) = f 



cPl 



T 
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dl 



2l + n 2 



xF(x + l + /p^ 



(2n-l+ 2 
k T + l T |) - 



n+l T ) 2 



Q{k 2 



T 



-F(x, 



(7) 



with a s = N c a s /ir, N c = 3 being the number of colors. We shall show that 
Eq. ([?D reduces to the conventional BFKL equation, the modified BFKL 
equation with the Q dependence, the modified BFKL equation with unitarity, 
and the new evolution equation with both the modifications by removing the 
applied approximations gradually. 

To derive the conventional BFKL equation, we simply assume the strong 
rapidity ordering, namely, approximate F(x + l + /p + ) by its dominant value 
F(x). Performing the integration over l + to infinity, and substituting F into 
Eq. (H), we arrive at 



dF(x, kr) 
dln(l/x) 



d 2 l T 

7rl T 



F{x, \k T + l T \) - 9{k 2 T - l T )F{x, k T 



(8) 



which is the BFKL equation adopted in [[|. 

However, the vanishing of F(x + l + /p + ) at large momentum fraction 
constrains / + to go to infinity. To render the approximation F(x + l + /p + ) ~ 
F{x) more reasonable, we truncate l + at the scale Q/V2 for the real gluon 
emission in Eq. (|7|), and obtain the Q-dependent BFKL equation, 



dF(x, kr) 
d\n(l/x) 



d 2 l 



a. 



-0£. 



T 



F(x, |k T + l r |) - 9(k T - l T )F(x, k T ) 
1t|) 



d 2 l T F(x, \k T 



7l 



P T + Q 2 



(9) 
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for the choice n = (1, 1,0). The extra term compared to the conventional 
BFKL equation comes from the upper bound of / + . It is trivial to find 
that Eq.(0) approaches Eq. (g) as Q — > oo. This correction, being negative, 
moderates the BFKL rise at low Q. Equation (§) has been solved in ||, and 
its predictions of the structure function F 2 match the data for various x and 

Q 2 . 

As stated in the Introduction, the conventional BFKL equation gives a 
power-law rise for the gluon distribution function at x — > 0, which violates 
unitarity. It has been pointed out that the assumption of the strong rapidity 
ordering is the cause for the unitarity violation For most values of 

l + , F(x + l + /p + ) is much smaller than F(x). Hence, replacing the former 
by the latter in the whole integration range of / + overestimates the real 
gluon contribution, which is responsible for the rise. To derive an evolution 
equation with unitarity, we employ Eq. (0) directly without applying the 
approximation. Reexpress F(x + l + /p + , \k T + l T \) as 



F(x, |k T + 1 T |) + F(x + l + /p + , |k T + 1 T |) - F(x, |k T + 1 T |) 



(10) 



in Eq. (|7|), where the first term, combined with the virtual gluon contribution, 
leads to the conventional BFKL equation, and the terms in the brackets are 
the correction from relaxing the strong rapidity ordering. We then derive 



dF(x, kx) 
d\n(l/x) 



a. 



+a 



d 2 l T 



F{x, |k T + Lr |) - 0(4 - %)F(x, k T ) 
4/+ 



(f 2 l - r x - 



dr 



X 



7T JO (2/+ 2 + /2)2 

F{x + l + / P + , |k T + 1 T |) - F(x, |k T + 1 T | 



:n) 



The above equation has been studied in |TD], and the BFKL power-law rise 
was found to be moderated into a logarithmic rise at small x. 



3. The new equation 

The Q-dependent BFKL equation @, though phenomenologically suc- 
cessful, still gives a power-law rise of the gluon distribution function. Hinted 
by Eq. ([IT]), we attempt to combine the above two modifications by switching 
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the upper bound of l + from oo to Q/V% in Eq. (|7|), and obtain 
dF(x, k T ) 



d\n(l/x) 



O.s 



rd 2 l T 


[1 


1 7T 





d v - 9^9^ ,2\2 F ( x + l kr + 1t D 



(y 2 Q 2 + i 2 T ) 



f fvrp Jjrp j 



Fix, k T ) 



(12) 



where the variable change l + = y^/xp + has been adopted, and the upper 
bound of y is determined by the kinematic relation Q « v / 2xp + . Equation 
( p!2|) is the new small-x evolution equation we shall investigate in more details 
below. 

In order to simplify the analysis, the 9 function for the virtual gluon emis- 
sion is replaced by 9(Qq — Ij) @, where the parameter Q can be determined 
from data fitting. This simplification is reasonable, because the virtual gluon 
contribution only plays the role of a soft regulator for the real gluon emission, 
and setting the cutoff of It to Qq serves the same purpose. We then Fourier 
transform Eq. (|l|) into the b space conjugate to kx with Eq. ([[(]) inserted, 
deriving 



dF(x, b) 
dWl/x) 



S{b,Q)F{x,b) 



+2a s {l/b)Qb / dyK 1 {yQb)[F{x^y^b)-F{x,b)] , 
Jo 



(13) 



where 



S(b, Q) = -2a s {l/b) lHQ b) + 7 - In 2 + K (Qb)} (14) 

comes from the combination of the first term in Eq. ( |10|) and the virtual 
gluon emission term. Note that S is the evolution kernel for Eq. ([|) in b 
space. K and Ki are the Bessel functions, and 7 the Euler constant. The 
argument of a s has been set to the natural scale 1/6. 

An initial condition F(xq, b) = F^ixQ, b) must be assumed when solving 
Eq. (|i~3D , xq being the initial momentum fraction. For instance, a "flat" gluon 
distribution function HI 



F (0 \x, b) = 3N g (l - xf exp(-Q 2 b 2 /4) 



(15) 
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for x > xo, N g being a normalization constant, has been proposed. Therefore, 
the initial function F^°'(x + y\fx, b) should be substituted for F(x + y\fx, b) 
in Eq. ( |T5D as x + y\fx > x . 

Before solving Eq. (|13D , we extract the behavior of F analytically. Insert- 
ing a guess F oc x~ x into Eq. (0), A being a parameter, we obtain 



\ = S + 2a s Qb I' dyK^yQb) 
Jo 



x + y\fx s 



x 



(16) 



It can be numerically verified that a solution of A, < A < S, exists for 
x < xq- That is, F increases as a power of x, consistent with the results from 
the conventional and Q-dependent BFKL equations. While the correction 
term, i.e., the second term on the right-hand side of Eq. ([16]), diverges as 
x —>■ 0, and no solution of A is allowed, implying that F can not maintain 
the power-law rise at extremely small x. We then substitute another guess 
F\ x ^o oc ln(l/x) with a milder rise into Eq. (|T3|) . In this case the correction 
term, increasing as ln 12 (l/x), only slightly dominates over the first term 
SF oc ln(l/x) (the derivative term —xdF/dx oc 1 is negligible), and Eq. (|13|) 
holds approximately. At last, we assume -F| x ^o ^ const, as a test. It is 
easy to find that the first term becomes dominant, and the correction term 
and the derivative term vanish, i.e., no const. ^ exists. These observations 
indicate that F should increase as ln(l/a;) at most when x approaches zero. 
In conclusion, the new evolution equation predicts a rapid power-like rise of 
F for x < Xq and a milder logarithmic rise at x — > 0. 

The structure function F2 is written, in terms of the /^-factorization 
theorem (IT2 



as 



F 2 (x, Q 2 ) = f 1 f T ^H(x/£, k T , Q)F& k T ) 

Jx t JO 71 



(17) 



p c being the upper bound of kr which will be specified later. The hard 
scattering subamplitude H denotes the contribution from the quark box di- 
agrams, where both the incoming photon and gluon are off shell by — Q 2 and 
— k\, respectively. H has been computed in 0, and its expression is given 
by 



H{z,k T ,Q) 



E 



2 a s 

e ^ z< 



z 2 + 



1-z) 



k 2 k 4 
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1 1 + Jl - 4z 2 k 2 /Q 2 ) 

x , In V . T/V -2 , (18) 

1 -4z 2 k 2 /Q 2 l-Jl- 4z 2 k 2 /Q 2 J 



with e g the electric charge of the quark g. To require a meaningful if, the 
upper bound of kx in Eq. flTT| ) is set to 

Pc = min (q, ^-Q^j . (19) 

We choose x = 0.1 and Q = 0.4 GeV, and evaluate Eq. ( |TTD for Q 2 = 15 
GeV 2 . Fitting the results to the corresponding data 0], the normalization 
constant is determined to be N g = 1.207. N g for other Q 2 are then extracted 
by requiring that the gluon density, defined by 

xg(x,Q 2 ) = f Q ^F(x,k T ), (20) 

JO 7T 

has a fixed normalization Jq xgdx. F 2 for Q 2 = 8.5, 12, and 20 GeV 2 are 
computed, and the results along with the data [0 are displayed in Fig. 1. 
For comparision, we also present the results from the conventional and Q- 
dependent BFKL equations 0. It is found that the shape of the curves 
from the conventional BFKL equation is almost independent of Q, and thus 
the match with the data is not very satisfactory. The curves from the Q- 
dependent BFKL equation, exhibiting smaller slopes for lower Q, match the 
data. However, they increase rapidly at x — ► 0, and violate the unitarity 
bound. The predictions from the new evolution equation also agree with 
the data well. The curves have a steeper rise at a larger Q, which is the 
consequence of the cutoff at the scale of order Q. Because of F(£, kx) < 
ln(l/£) at £ — >• 0, F 2 rises as ln 2 (l/x) at most as indicated by the £ integration 
in Eq. flTTp, and thus satisfies the unitarity bound. It is observed that the 
ascent of the curves does not speed up as obtained from the conventional 
and Q-dependent BFKL equations at x < 10~ 4 . This is the consequence of 
relaxing the strong rapidity ordering. 



4. Conclusion 

In this letter we have combined the modifications to the conventional 
BFKL equation obtained in our previous studies, which are based on the 
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resummation formalism. The resultant small-x evolution equation possesses 
the intrinsic Q dependence and satisfies the unitarity requirement. The pre- 
dictions of the structure function F 2 are in agreement with the HERA data 
for various x and Q 2 . We have also observed that the power-law rise turns 
into a logarithmic rise as x < 10~ 4 . 

We emphasize that the cutoff of order Q for the longitudinal momentum 
in the real gluon emission is phenomenologically motivated in this work. In 
the derivation of a new unified evolution equation for both intermediate and 
small x, we shall demonstrate that the scale Q can be introduced in a rigorous 
way. This subject will be published elsewhere JL3|. 

This work was supported by the National Science Council of R.O.C. under 
the Grant No. NSC-87-2112-M-006-018. 
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Figure Caption 

FIG. 1. The dependence of F 2 on x for Q 2 = 8.5, 12, 15, and 20 GeV 2 
from the conventional BFKL equation (dotted lines), from the Q-dependent 
BFKL equation, and from the new evolution equation (solid lines). 
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